Entanglement of transverse modes in a pendular cavity 



Stefano Mancini and Alessandra Gatti 
INFM, Dipartimento di Fisica, Universita di Milano, Via Celoria 16, 1-20133 Milano, Italy 

(February 1, 2008) 

We study the phenomena that arise in the transverse structure of electromagnetic field impinging 
on a linear Fabry-Perot cavity with an oscillating end mirror. We find quantum correlations among 
transverse modes which can be considered as a signature of their entanglement. 
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■ I. INTRODUCTION 

o ' 
O ' 

I It is now well assessed that an empty optical cavity with a moving mirror in its steady state may mimic a Kerr 
^ . medium if illuminated with coherent light JfJJ^] . When the mirror is free to oscillate, the radiation pressure induces 
a coupling between its position and the intensity of the light beam, modifying the optical path in an intensity 
dependent way. A wide range of applications of this effect have been recentely developed (3|. In particular the 
^\ ' system, showing a typical bistable behaviour, can be used as a quantum noise eater device, because the output light 
t-H . is significantly squeezed dH , or to generate highly nonclassical states for both radiation and mirror HQ , due to the 
nonlinear character of the interaction. The theoretical treatments of this model were fully quantum but at most one 
dimensional in space because of the plane wave approximation, which ensures that the electric field is uniform over 
the transverse plane. This means that the investigations dealt only with temporal-frequency aspects, neglecting all 
features related to space. On the other hand, recent years have seen an increasing interest towards the spatial aspects 
|^| . Hence, the aim of this paper is to study the spatial phenomena that arise in the transverse structure of the 
light field impinging on a linear Fabry-Perot cavity with an oscillating end mirror. We shall show the possibility to 
correlate transverse modes, other than to obtain squeezing effects. Such correlations have different nature, and could 
reveal entanglement ]9|,[l0|]. Moreover, we will analyse the differences and analogies with a Kerr nonlinear system, and 
we will investigate the influence of temperature on the transverse structures for such optomechanical system. 
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' II. THE MODEL 

We consider a linear Fabry-Perot empty cavity with one fixed mirror, partially transmitting, and one perfectly 
reflecting end mirror. The completely reflecting mirror, having a mass m, can move, back and forth along the cavity 
^ , axes (say z), undergoing harmonic oscillations at frequency uj m . The amplitude of such oscillations is however, much 
less than the equilibrium cavity length L. The cavity resonances are calculated in absence of the impinging light. The 
characteristic cavity frequencies are assumed many orders of magnitude greater than u> m to ensure that the number 
of photons generated by the Casimir effect is completely negligible. We also assume that the cavity round trip time 
is much shorter than the mirror's period of oscillation and the Doppler frequency shift of the photons jllj on the 
moving mirror, is completely negligible. 



A. The Single Mode Model 

Let us recall the model describing the system in the approximation where all the spatial effects (e.g. diffraction) 
are neglected and the cavity is assumed to operate in a single longitudinal mode of frequency u)q. That cavity mode 
is coherently driven by the action of an external field £e~ lu)at , having a frequency u s and a complex amplitude £. 
By indicating with A, A^ the annihiliation and creation operators of the cavity mode, the system Hamiltonian, in a 
frame rotating at the frequency w s , reads Q 

H = hAA^A + huj m (— + — )- hgA^AX + ih {£A ] - £* A) , (I) 

where A = loq — oj s is the cavity detuning; X and Px are the dimensionless position and momentum operators of 
the mirror, obeying the commutation relation [A, Px] = i. The interaction part (third term) of Eq.(IH) accounts for 
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the effect of radiation pressure force which causes the instantaneous displacement of the mirror, and the coupling 
constant is given by || 



L V rnuj m 

In writing down the equations describing the dynamics of the system, we must take into account the damping of 
the movable mirror due to the coupling with a thermal bath in equilibrium at temperature T, and the cavity losses, 
due to the coupling of the internal mode with all the external modes of radiation through the fixed (transmitting) 
mirror. Hence, we have the following quantum Langevin equations 

d t A{t) = -iAA(t) + igX{t)A{t) + E- lc A(t) + y/2^ c A in (t) , (3) 
9tX{t)=u m P x {t), (4) 
d t Px(t) = -uj m X(t)+gA^(t)A(t) - lm Px{t) + y/l^e in {t) , (5) 

where j c , 7 m are the decay rates of the cavity mode and of the mirror momentum, respectively. The noise operators 
(labeled with the subscript "in") have zero expectation value, and obey the following correlations: 

(A in (t)A in (t>)) = , (AUt)AUt')) = S(t - t') , (6) 



(e in (t)e in (t'))=N T S(t-t'), (7) 

where Nt is the number of thermal excitations of the mirror Nt = fcsT/?iw m , with ks being the Boltzmann constant. 
The form of Eq.(^) with correlations (Q) corresponds to the quantum optical master equation JL3]. Instead, the form 
of Eqs.(|4j), (||), with correlations (^), corresponds to the standard quantum Brownian master equation [^2). The latter 
is only valid in the limit Nt 3> 1, while a more careful analysis is required in the opposite limit [pi. 



B. The Spatially Multimode Model 

Motivated by the recent progress in the study of transverse quantum effects in optical systems |§ we wish to extend 
the previous model to the case of spatially multimode fields. Let us again assume the validity of the single longitudinal 
mode approximation, which, toghcthcr with the mean field approximation allows us to neglect the dependence of the 
field over the variable z. However, we now allow the radiation field to depend on the transverse vector x = (x,y), 
which is the position vector in the plane orthogonal to the direction z of propagation of fields. 

Then, Eqs.(||), (^), (||) after elimination of Px, must be rewritten as 

8 t A(x, t) = -»AA(x, t) + igX(t)A(x, t) + £(x) - 7c A(x, t) + i lc P D ^\M*, t) + VWc A„(x, t) , (8) 
3 2 t X{t) = -u 2 m X{t)+gLu m j dxAt(x,t)^(x,t) - j m d t X{t) + (*), (9) 

where the transverse Laplacian = d 2 /dx 2 + d 2 /dy 2 has been introduced to describe the diffraction in the paraxial 
approximation [Q. The coefficient l 2 D has the dimension of an area, and is given by £ 2 D = c 2 /2ui s "f c . It represents 
the typical length scale for spatial structures emerging in an optical resonators. In the following we shall rescale all 
the spatial lengths over Id- 

For the sake of simplicity we considered only the case of a uniform mirror motion over the transverse plane, even 
if there is the possibility to take into account acoustic modes of the resonators but this is planned for a future 
work. 

In order to avoid difficulties arising from the continuum of transverse modes, we consider in the transverse plane 
(x,y) a square of side t, and we assume periodic boundary conditions for the fields. A complete set of transverse 
modes (corresponding to a single longitudinal resonance) is then given by 

1 2tt 
/n(x) = - exp («k n • x) , with k n = — n, n=(n x ,n y ), (10) 

where (n x , n y ) is a couple of integer numbers, n x ,n y = 0, ±1, ±2 . . .. Then, the fields can be expanded in the following 

way 
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A(x, t) = £ /n(x)a n (t) =jJ2 e * n ' Xa n(0 , 
n n 

f(x)=^/ n (x)e n = i^e*" x en . 

1 1 n 

The Hamiltonian (Q) of the single mode model can thus be generalized as 

H = huj m ( + — J + ^ [hAna^an] - hg ^ [ai,a n l] + [ a n e n - a n e* ] , 

^ ' n n n 

where we introduced the mode detuning A n = uj n — lo s , with 

the frequency of the transverse mode n. 

We are now able to derive from Eq.(^) the Langevin equation for each transverse mode 

d t a n (t) = -iA n a n (t) + igX{t)a n {t) + e n (t) - a n (t) + V2a™(t) , 
dfX(t) = -ujf n X(t) +5w m ^a| n (t)a m (t) - ■y m d t X(t) + Lu mV ^e in (t) , 



(11) 
(12) 

(13) 
(14) 



(15) 
(16) 



where we introduced the following scalings 
J c t->t, U} m /lc^ u m , Im/lc^lm, dhc^-g, £/lc->£, A n /7 c ->A n , Ai n /^/7c — > A in , ejn/^/Tc — ► Ein ■ 

Equations (|l5|), ( |l6| ) show that each mode interacts with the mirror, which in turn becomes an intermediary between 
the various modes, redistributing the quantum information among them. Depending on such a process we expect that 
the output light will have different characteristics from the input one. 

III. THE STEADY STATE 



We are actually interested in the classical steady state regime and in small fluctuations around this steady state. 
To this end, as it is usually in the semiclassical treatment of quantum noise, we set a n = a n + <5a n , X = x + 5X, 
where the c- numbers a n , x are solutions of the classical steady state equations 



= — iA n a n + igxa n + e n — a n . 
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(17) 
(18) 



By eliminating the variable x from Eqs.(17), (|18|), it is possible to get an infinite set of coupled cubic equations 
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a n < 1 + i 



Ei 



where we have scaled the variables as follows 

2 \ 1/2 
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(19) 



(20) 



In the case of only one spatial mode, the above set of equations reduces to only one cubic equation which shows a 
typical bistable behavior ||. 

To proceede on we consider, without loss of generality, a gaussian and real pump 



V 7T Wp 

where w p idicates the pump waist (we will always use i ^ w p ). The total input power is given by 



(21) 
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while the pump components on the spatial modes are 



J d X \£( X )\\ 



t — - exp 



Now, by defining the intracavity power P as 

P = J dx|a(x)| 2 =^|a 



2 

n| ; 



(22) 



(23) 



(24) 



from Eq. (091) we can obtain 



l«n| 2 = 



1 + (P-A n )2 ' 



(25) 



and summing over the index n 



Pi, 



He 



27r(^) 2 exp[-(^)kl 



1 + [A„ - py 



(26) 



This equation gives the functional relation between the incident and intracavity intensity. The bistable behavior for 
the total intensity is shown in Fig.l. Depending on the slope of the curve, we have stable and unstable branches. 
Once one has choosen the working point along this curve, the input and intracavity powers are well defined, hence it 
is possible to calculate each single steady state components of the intracavity field as 



<Pw 



^ ^ exp 



1 + i [A n - P] 



(27) 



By considering also the field reflected at the transmitting mirror we have the input-output relation JL6| , in terms 
of scaled variables 



cw(x) = 2a (x) - a in (x) 



(28) 



where £(x) = a; n (x). It is possible to see, in Fig. 2, that the Gaussian profile of the beam is maintained from input 
to outgoing fields. 

Finally, the stationary displacement of the mirror due to the radiation pressure results x = P/g. 



IV. DYNAMICS OF SMALL FLUCTUATIONS 



The evolution equations for small fluctuations coming from the linearization of the Eqs.(|l5|), ( |l6| ) around a stable 
steady state are 



d t Sa n {t) = -iA n 5a n (£) + iy/iJZSX(i)a n + iP 5a n (t) - 5a n (t) + V25a™(t) , 

dl 5X(t) = -w^SX^) + w m ^/u^2_j [ a ™^ a ln(*) + o* m 8a m {t)\ - j m d t SX(t) + £m(t) 



(29) 
(30) 



We immediately recognize that in the case of fiat pump, i.e. a n = ao <5n,o, the coupling with the mirror occours only 
for one mode (the fundamental), hence no transverse effects can arise. This is one of the peculiar differences with 
respect to the Kerr-like model, and it is due to the fact that only one mode practically survive at steady state. 
In the frequency domain Eqs.(^), ( (30| ) become 

(31) 
(32) 



in8a a (n) = -iA n (5a n (fi) + iy/uJ^SX^a^ + iP <5a„(ft) - <5a n (ft) + V2Sa™(n) , 



k 
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By eliminating the mirror variables one gets an infinite set of linear equations 



[iCl + iA n -iP + l] Sa a (Cl) - ia nX {Cl) [a^Sa^-Cl) + a* m Sa m (Ct)] = y/2 8a£(il) + ia n x(P)J— , (33) 



where we have introduced the mirror response function 

Other typical differences between this model and a Kerr nonlinear system are due to the dynamics of the moving 
mirror characterized by a frequency dependent susceptibility, and to the presence of a thermal noise as can be seen 
from Eqs.^3) and (^). 

In Eq. (|33|) we may see that the various radiation modes can become correlated and some spatial effects should 
appear. The latter should also depend on temperature. 



V. SOLUTIONS 



In order to find the solutions of Eqs.(B3J), let us write them in a matricial form. In doing so, we introduce a truncation 
in the number of modes effectively achieved, namely we let n Xl n y = — n, . . . ,n. The truncation is reasonable given 
the fact that the pump field supports a finite number of modes. 

Then we introduce the following vectors 



V 



A 



5hl -(-CI) 



5a-n t n(Cl) 
5a" 1 (-CI) 

6an y n(Cl) 

*4n(-n) 



They all have 2(2n + l) 2 components. 

Now, the system of Eqs.(p3[), can be written as 



lot— n,— n 



—ia% — 



-ia%-- 



V 



^ — n. — n 
^ — n. — n 



^ n. . — r 



(35) 



[M + (iCl + 1)21 • V = V2 V in + . -^ x (tyAi in (Cl) 



(36) 



where the vector V m is defined analogously to Eq.pf]), but for the input operators. Instead, X is the 2(2n + l) 2 x 
2(2n + l) 2 identity matrix, while M is the 2(2n+ l^x 2(2n + l) 2 matrix given by 



M k>l = t(-) k [P - V k ] 5 kil + i(-) 1 X (Cl) A k A-(-i)' 
By inverting the relation ( |36| ) we can get the formal solution of the system as 

V2V in + <f^x(n)Ai in (Cl) , 



V = T ■ 



F=[M + (iCl + l)!}- 1 



(37) 



(38) 



However we are looking for the solution of the outgoing modes. Then, the input-output relation [ fig ] can be written 
in vector form as 
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and if we combine it with Eq.(p8[), we get 



where 



V out = V2V- V" 



V out =B- V m +Ue h 



B = 2 T -1. 



and 



(39) 
(40) 
(41) 

(42) 



From the above matrix relations one can extract the expression for the various components 8a^ n (f2) , 55°"*^ (f2) 
in terms of the input noise operators. 

Finally, it is worth noting that the transformation among input and output fields could not preserve the commutation 
relations. This can be understand by observing Eq.(|37j) where Ke{Mk,k} yields and additional damping term 7 a dd oc 
Im{x(^)}- Neverthless, if 7 a dd "C 1 the commutation relations are preserved. That happen for istance in the case of 

n -> o. 



VI. THE OUTPUT CORRELATIONS 



Using Eq. (^) , and the fact that all the elements of the vector V ln are uncorrelated except those of the form 
we are able to calculate the correlations of the output modes from the components of Eq.(Eol). They result: 



(27T+1) 2 



B{n}+2,2k-l B{n} + l,2k + N T U{n}+2 ^{n} + l 



fc=l 



5(n + n') 



and 



(2n+l) 2 

^ S{n}+l,2fc-l 6{n}+l,2fe + N T ^{n} + l ^{n} + l 



k=l 



6{Cl + il') , 



(43) 



(44) 



(45) 



where {n} = 2(2n + l){n y + n) + 2(n x + n). 

We now introduce the linearized output intensity operator 

SI° ut (n) = a n ut *S5, n ut {Q) +< ut <5a° utt (-f7) . 

Then, the spectrum of the output intensity correlations among the various modes, namely 



S2%(n) = J dn> 



(46) 



(47) 



1 1 i in 



can be easily calculated by using Eqs.([44|), ((45|). It should be compared with that of a coherent state, i.e. ^"^(O) = 

Furthermore, by considering the total output intensity SI otLt (£l) = X)n^n Ut (^)i it i s possible to define the output 
intensity spectrum 



s out {n) = J dn 



, (s r ut (n)si out (op) 



(48) 



which, again, can be easily calculated by using Eqs.(Ej), (f45[). It should be compared with that of a coehrent state, 
i.e. 5 0Ut (0) = 1. 
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VII. RESULTS AND CONCLUSIONS 



To study numerically the system, we have considered a pump defined by (19 x 19) modes. 

In Fig. 3 we show the intensity (auto)correlation spectrum for each mode, S^^O = 0.1). The squeezing is found in 
several modes around the fundamental, but it tends to disappear when the temperature increases. The total intensity 
squeezing 5 out (fi = 0.1) has been calculated for the three cases of Fig. 3 obtaining the following values: 0.11, 0.62, 
and 1.14. The order of magnitude is similar to that obtained in Refs. H,||, however, in this case the total squeezing 
is distributed among the various mode. 

In Fig. 4 we plot the intensity correlation spectrum of the fundamental mode with each other 5o u „(f2 = 0.1). At 
low temperature we have a negative correlation of the fundamental mode with the neighbourhoods (the central peak 
instead represents the autocorrelation of the fundamental mode, i.e. its intensity squeezing). The negative correlations, 
i.e. anticorrelations, can be understood by considering the interaction with the mirror as redistributing the photons 
from the fundamental (pumped) mode to its neighbourhoods. On the other hand these anticorrelations could be 
considered as a signature of entanglement. In fact, thought the notion of entanglement is not clear for open system, 
we can see from Eqs.(^), ( Kj] ) that the correlations among various modes arise as consequence of both vacuum noise 
and thermal noise. However, the latter, which is of classsical origin, only leads to positive correlations. Instead, the 
vacuum noise, purely quantum, can give positive as well as negative correlations. Therefore, anticorrelations can 
be considered as a signature of purely quantum correlations, hence, entanglement. Nevertheless, the two types of 
correlations become competing as can be seen from top to bottom of Fig.[|. As matter of fact, by increasing the 
temperature the negative correlations tend to disappear, and the shape tends to assume the form of <5o,o- The residual 
(positive) correlations may be attributed to the thermal noise and could be used to study the Brownian motion of 
the mirror jPjJ. It is to remark, anyway, that the quantum effects are quite robust to the thermal noise, providing 
to have an high quality factor for the moving mirror uj m /lm- This can be easily understood by noticing the factor 
multiplying the thermal noise term in the set of equations (|33|). 

Thought we limited our analysis to the intensity correlations, we might argue the existence of more fundamental 
correlations, like EPR correlations |fl8| . 

In conclusion, we have studied the phenomena related to the finite extend of a light beam in optomechanical 
coupling. In some sense this work can be considered complementary to Ref. E^), where instead several vibrational 
modes of the mirror were coupled to only one light mode. Moreover, the system can give the possibility of multimode 
entanglement which is one of the most striking aspect of quantum mechanics PJTfj| . Finally, the developed theory 
could be also used in gravitational interferometry ]l9| where such transverse effects could increase the measurement 
sensitivity. 
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FIG. 1. The bistability curve is represented for several values of detuning. The chosen pump waist is w p — 2. 



Re [a m] 




-n . 4 1 



Re [a out] 

FIG. 2. The real and imaginary part of the input and ouput stationary fields are represented as a function of r = (x 2 + y 2 ) 1 ^ 2 
for w p = 2, P in = 2.89, P = 1.06, A = 2. 
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FIG. 3. The spectrum S°"*{Q = 0.1) is plotted vs n for several values of N T (from top to bottom iV T = 10 4 , 10 5 , 10 6 ). 
Values of other parameters are: w p = 2, Pi n = 2.89, P = 1.06, Ao = 2, mechanical quality factor u) m /"j m = 10 6 . 
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FIG. 4. The spectrum Sg^i(fi = 0.1) is plotted vs n for several values of N T (from top to bottom N T = 10 4 , 10 5 , 10 6 ). 
Values of other parameters as in Fig. 3. In reality, the central peaks of Fig. 4 go out of the actual scale, and their values coincide 
with those of Fig. 3. 
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